A SuperMassive Black Hole Fundamental Plane for Ellipticals 

Sudhanshu Barway 1 and Ajit Kembhavi 2 

Inter University center for Astronomy and Astrophysics, Post Bag 4, Ganeshkhind 

Pune 411 007, India 

ABSTRACT 

We obtain the coefficients of a new fundamental plane for supermassive black 
holes at the centers of elliptical galaxies, involving measured central black hole 
mass and photometric parameters which define the light distribution. The galax- 
ies are tightly distributed around this mass fundamental plane, with improvement 
in the rms residual over those obtained from the Mbh — o~ and Mbh — L rela- 
tions. This implies a strong multidimensional link between the central massive 
black hole formation and global photometric properties of elliptical galaxies and 
provides an improved estimate of black hole mass from galaxy data. 

Subject headings: black hole physics - galaxies: elliptical and lenticular, cD - 
galaxies: nuclei - galaxies: evolution - galaxies: fundamental parameters - galax- 
ies: kinematics and dynamics 

1. Introduction 

The existence of massive black holes (hereafter BH) at the center of nearby inactive 
galaxies, as well as in the nuclei of active galaxies and in quasars, is well established. Obser- 
vations based on high resolution data and reverberation mapping are now available which 
allow measurement of the masses of BH using different techniques (Ferrarese & Ford 2005; 
Metzroth et al. 2006; Shapiro et al. 2006). Kormendy & Richstone (1995) showed that the 
measured BH mass M B h is correlated with the bulge luminosity L and bulge mass M bu i ge 
with rms scatter ~ 0.5 dex in log Mbh (see also Magorrian et al. 1998). A tight correlation 
between Mbh and the central velocity dispersion a of the host galaxy with smaller rms scat- 
ter of ~ 0.34 dex in log Mbh was reported by Ferrarese & Merritt (2000) and Gebhardt et 
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al. (2000); however, the published estimates of slope in M B h- o relation span a wide range 
(3.75-5.30, see Tremaine et al. 2002). The small scatter of the Mbh- cr relation suggests 
that the bulge dynamics (or mass), rather than the luminosity, is responsible for the tight 
correlation. 

It is believed that massive black holes play an important role in the formation and 
evolution of galaxies, and the growth of the BH and bulges must be linked to the same 
physical processes; this results in BH masses that are related to the properties of host 
galaxies (Silk & Rees 1998; Haehnelt & Kauffmann 2000; Adams et al. 2001; Merritt & 
Poon 2004; Sazonov et al. 2005). Graham et al. (2001) and Marconi & Hunt (2003) have 
shown that when bulge parameters are measured with sufficient accuracy using the technique 
of bulge-disk decomposition, the resulting scatter in the M B n-L relation is comparable to 
that in the M Bii -a relation (see also Graham 2007). Marconi & Hunt (2003) also suggested 
that a combination a and bulge effective radius r c should be used to derive the correlations 
between Mbh and other bulge properties. Recently, Lauer et al. (2006) have suggested that 
the bulge luminosity may be a better indicator of BH mass than the bulge velocity dispersion 
at the high mass end for brightest cluster galaxies. However, in spite of all these attempts, 
our understanding of how the photometric properties of galaxies and their central BHs are 
linked in the process of formation of galaxies remains unclear. 

In this Letter, we show that log Mbh, logr e , and (fib{< r c)), which is the mean bulge 
surface brightness in magnitude within r e , are tightly correlated for nearby elliptical galaxies 
having measured central BH masses. The scatter around the best fit plane is significantly less 
than the scatter in various two-dimensional relations. It is also less than the scatter obtained 
if BH masses are estimated from the photometric parameters of galaxies using the standard 
fundamental plane for ellipticals and the M B h-o" relation. In §2 we provide details about the 
samples of galaxies used in the analysis. We present the results in §3, a discussion in §4 and 
in §5 a summary of the work. Throughout this Letter, we use H = 70 km sec -1 Mpc -1 , and 
express r e in kiloparsec, a in units of kms -1 , and mass and luminosity in Solar units. 

2. The Data 

To obtain the photometric scaling relation we have considered a sample of 20 galaxies 
classified as elliptical in the Ferrarese & Ford (2005) galaxy list with measured black hole 
masses. In Table 1 we report the relevant data for this sample. To compare the estimates 
of central black hole masses obtained from our planar relation and the Mbh- ct and Mbh-L 
relations, we consider a sample of 22 elliptical galaxies from the Coma cluster. This sample 
was observed by Jorgensen et al. (1992) in the Johnson B band; a description of the data 
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and the global parameters obtained from the images can be had from the reference. 

3. A New Fundamental Plane for Nearby Ellipticals 

The Mbh- ct arid Mbh- L relations offer two ways to estimate the BH mass from other 
galaxy properties, and have been applied to AGN (McLure & Dunlop 2002), BL Lac objects 
(Falomo et al. 2002), low-redshift radio galaxies (Bettoni et al. 2003) and to bright cluster 
galaxies (Lauer et al. 2006; Batcherdor et al. 2006). We have revisited the Mbh- ct relation 
and Mbh- L relation by applying a bisector linear regression fit (Akritas & Bershady 1996) 
to the data given in Table 1 for the sample of nearby elliptical galaxies with measured BH 
masses. The two best fit relations are: 

log M BH = (4.53 ± 0.49) loga- (2.24 ± 1.17) (1) 

log M BH = -(0.56 ± 0M)L B - (3.10 ± 1.51) (2) 

The rms scatter around the best fit lines above is 0.34 dex and 0.42 dex respectively, along 
the log Mbh axis. Both the relations are in good agreement with those in Bettoni et al. (2003) 
and reference therein, but the relations are different from those of Ferrarese & Merritt (2000) 
and Gebhardt et al. (2000), as we have used a sample of nearby ellipticals only. It is possible 
that some of the scatter seen in Mbh — & relation and Mbh — L relation is caused by the 
effect of a third parameter. This is supported by the strong correlation that we find between 
log Mbh and logr e , with a correlation coefficient r = 0.89, which is significant at the 99.99 
% confidence level for 19 objects; Marconi and Hunt (2003) have obtained a similar result. 

Our aim is to derive a planar relation involving the BH mass and the basic photometric 
parameters r e and (/Ub(< r e)); this can be used to estimate the black hole mass when it is not 
known from measurement, without reference to a spectroscopically measured quantity like 
the central velocity dispersion. We find that the least scatter around the best-fit plane in the 
space of the three parameters is obtained by expressing it in the form log r e = a log M B h + 
b (/Ub(< r e )) + constant. We minimize the sum of the absolute residuals perpendicular to the 
plane, excluding one galaxy NGC4742, which is an outlier we have identified in Figure [TJ 
The equation of the best-fit mass fundamental plane is 

logr e = (0.32 ±0.06) log M BH + (0.31 ±0.06) (fi h (< r e )) 

-8.69 ±1.58 (3) 

The uncertainties on the mass FP coefficients were determined using a bootstrap method. 
An edge-on view along log r e of the plane is shown in Figure [T](a). The rms scatter in the 
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direction of logr is 0.061 dex. Figure [T](b) shows another edge-on view of mass FP in the 
direction of log Mbh, with rms scatter in that direction of 0.19 dex, which is significantly 
less than the scatter in the M^-a relation (Gebhardt et al. 2000; Ferrarese & Merritt 2000). 
The outlier NGC4742 is 6.32 x (rms scatter) from the plane along the logr e axis. We have 
also obtained the equation of the best fit plane including this outlier. The rms scatter then 
increases to 0.078 dex in log r e and 0.25 dex along log Mbh axis respectively. Therefore, even 
with the outlier included we have less scatter than in the log M B h — log o and log M B h — log L 
fits. 

If we exclude from the fit the four galaxies NGC 821, NGC 2778, NGC 4649 and NGC 7052 
for which the BH sphere of influence is not resolved, and the outlier from the fit, the rms scat- 
ter in log M B h around the best-fit plane obtained using the remaining 14 galaxies decreases 
to 0.17 dex. 

For nearby ellipticals we have derived the standard fundamental plane relation, using 
the same technique as in the case of the mass fundamental plane and again excluding the 
outlying data point NGC 4742. The equation of the best-fit FP is 

logr e = (1.34 ±0.22) log a + (0.30 ±0.05) (/x b (< r e )) 

-8.93 ±0.74 (4) 

The rms scatter is 0.068 dex in log r e . The FP coefficients and rms scatter around the fit 
are in agreement with those available in the literature (Jorgensen et al. 2006). 



4. Discussion 

As suggested by Ferrarrese & Ford (2005), given the photometric parameters of an 
elliptical galaxy, the central velocity dispersion a can be derived using the FP relation given 
in Equation HI if it is not directly observed, and then the Mbh- 0" relation in Equation [1] can 
be used to estimate the BH mass. However, the error in the estimated BH mass will then 
be the cumulative error of these two relations, thus increasing the uncertainty in the mass 
estimate. Another disadvantage of this approach is that the slope in the Mbh — cr relation 
spans the range 3.75-5.3, leading to further uncertainty in the estimate of the mass. The 
mass FP provides an improvement over this two step procedure, and also helps to constrain 
the slope of the M B h- & relation, as described below. 

We consider a two-dimensional relation of the form log Mbh = a logo" ± /?, where a 
and /3 are constants to be determined. Introducing this into Equation H] for the fundamental 
plane, we get a plane in the space of log Mbh, logr e and (/it>(< r c )), with the direction of 
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the normal to the plane dependent on the value of a. In Figure [2] we have plotted, as a 
solid line, the angle between this normal, and the normal to the mass FP in Equation [31 for 
a range of values of a. The filled circles on the curve indicate the angles corresponding to 
specific values of a found in the literature, obtained by various groups from their fits to the 
data (see Tremaine et al. 2002). It is seen from the figure that the angle between the two 
planes is minimum near a = 4.5, which should be the value to be used in the log Mbh — log cr 
relation to determine black hole mass from the central dispersion velocity. The best fit in 
Equation [T] corresponds to a = 4.53. It will be interesting to see how the minimum value of 
a depends on the morphological type of the host galaxy. 

We have used the mass FP to predict the black hole mass for a set of 22 elliptical 
galaxies from the Coma cluster, using photometric data from Jorgensen et al. (1992). We 
have also obtained the black hole mass for these galaxies using Equations [1] and [2j The masses 
obtained in these various ways are compared in Figure [3j It is seen from Figure [3]^a) that the 
agreement between M B H(mass FP) and M B n(a) is good; the points are distributed around 
a line with slope close to unity, with a correlation coefficient r = 0.93, which is significant 
at better than the 99.9 % level. A larger number of points will be needed for a better 
comparison and to examine any departures from linearity. The slope in the log Mbh — log cr 
relation in Equation [1] is close to the minimum value of a obtained from Figure [2j Using 
any other value of a will produce a less favorable comparison. We see from panel (b) that 
Mbh (mass FP) and Mbr(L) are distributed along a straight line with slope less than unity; 
for Mbh ~1O 8 ' 5 M , masses obtained from the FP would be systematically less than masses 
obtained from the log Mbh — log L relation. The dispersion of the points around the best 
fit line is greater in this case than in panel (a). We have for completeness compared in 
panel (c) black hole masses obtained from the log Mbh — log cr and log Mbh — log L relations 
respectively, and find a slope greater than unity and larger dispersion than in the other cases. 

The three dimensional mass FP has lesser rms deviation than in the earlier two dimen- 
sional relations while some reduction in residuals is expected when the number of parameters 
in the fit is increased from two to three, it appears that the process can not be taken any 
further. We have considered a four dimensional plane with the dispersion velocity a included 
in the fit along with the two photometric parameters. However, we find that the residuals 
from the three dimensional plane are not correlated with logcr, and the quality of a four 
dimensional fit involving log Mbh, logr e , (Hb(< r c)) and log a is poor. A three dimensional 
relation is therefore the best we can do with the available data. 

It will be interesting to obtain the mass FP for photometric data in the near-infrared 
bands, since stellar population metallicity effects are less important than in the optical region 
(Pahre et al. 1998). Another issue to examine is whether the bulges of galaxies of various 
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morphological types share a common mass FP. 



5. Summary 

We have shown that logr e , logMBH and (^b(< r e )) for nearby elliptical galaxies having 
measured central BH masses are tightly distributed about a plane with a rms scatter of 0.19 
dex along logM B H- The scatter decreases to 0.17 dex in logM B H when we use only those 
galaxies for which the BH sphere of influence is resolved. The mass FP provides a convenient 
way for estimating BH mass from photometric data alone. 
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Table 1. Basic parameters for elliptical galaxies with measured black hole mass. 



Object 


Type 


Distance 






a 


L B 


log r e 


(Wb(< r c )) 






(Mpc) 


(10 s M ) 


(km s _1 ) 


(mag) 


(kpc) 


(mag arcsec -2 ) 


NGC 221/M32 


-6.0 


0.80 


9 c:+0.5 
z -°-0.5 x 


10 6 


75zbl0 


-15.80±0.18 


-0.83 


18.69 


NGC 821 


-5.0 


24.1 


o 7 +2.4 
°-'-0.8 


10 7 


209±26 


-20.42±0.21 


0.72 


21.85 


NGC 2778 


-5.0 


22.9 


1.4t°J x 


10 7 


175±22 


-18.58±0.33 


0.26 


21.38 


NGC 3377 


-5.0 


11.2 


1.0+g"? x 


10 8 


145zbl7 


-19.18±0.13 


0.26 


20.76 


NGC 3379 


-5.0 


10.6 


1.01™ x 


10 8 


206±26 


-19.81±0.20 


0.26 


20.16 


NGC 3608 


-5.0 


22.9 


i-s±S:8 >< 


10 8 


182±27 


-20.07±0.17 


0.59 


21.41 


NGC 4261 


-5.0 


31.6 


5.2+1° x 


10 8 


315±38 


-21.23±0.20 


0.77 


21.25 


NGC 4291 


-5.0 


26.2 


3.1+0"! x 


10 8 


242±35 


-19.72±0.35 


0.27 


20.25 


NGC 4374/M84 


-5.0 


18.4 


i-ol3:8 x 


10 9 


296=b37 


-21.40±0.31 


0.68 


20.81 


NGC 4473 


-5.0 


15.7 


l.ll°i X 


10 8 


190±25 


-19.86±0.14 


0.28 


20.19 


NGC 4486/M87 


-4.0 


16.1 


3.4+;° x 


10 9 


375±45 


-21.71±0.16 


0.91 


21.60 


NGC 4564 


-5.0 


15.0 


5.6+8J x 


10 7 


162±20 


-18.94±0.18 


0.19 


20.64 


NGC 4697 


-5.0 


11.7 


1.71°:? x 


10 8 


177±10 


-20.20±0.18 


0.63 


21.41 


NGC 4649/M60 


-5.0 


16.8 


2-0+0;| x 


10 9 


385±43 


-21.30±0.16 


0.78 


21.10 


NGC 4742 


-5.0 


15.5 


1-4181 x 


10 7 


90=b05 


-19.03±0.10 


-0.06 


19.36 


NGC 5845 


-5.0 


25.9 


2.4l°+ x 


10 8 


234±36 


-18.92±0.25 


-0.30 


18.38 


NGC 7052 


-5.0 


71.4 


4.0+ 2 | x 


10 8 


266±34 


-21.43±0.38 


0.89 


22.01 


IC 1459 


-5.0 


29.2 


i.5+;° x 


10 9 


340±41 


-21.45±0.32 


0.73 


20.81 


NGC 6251 


-5.0 


107.0 


6.1 + |° x 


10 8 


290±39 


-21.95zb0.28 


1.31 




CygA 


-5.0 


240.0 


2.9±0-7 x 


10 9 


270=b90 


-20.03zb0.27 







Note. — Cols. 1 and 2 give the name and the morphological type from RC3; Col. 3 the distance, derived 
from Surface Brightness Fluctuations (SBF, Tonry et al. 2001); Cols. 4-6 provide the adopted values for 
the mass of black hole Mbh> velocity dispersion and absolute bulge luminosity L in B band (from Ferrarese 
& Ford 2005); Cols 7 and 8 give the effective radius r e (from Faber et al. 1989 and using the distance in 
Col. 3) and mean surface brightness within effective radius in B band (from Faber et al. 1989). 
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-1 1 2 5 6 7 8 9 10 11 

0.32±0.06 log M BH + 0.31±0.06 </i„(<r e )> - 8.69±1.58 3.17±0.48 log r, - 0.99±0.22 <M b (<r„)> + 27.54±4.36 



Fig. 1. — Edge-on views of the mass fundamental plane relations for nearby ellipticals: (a) 
along one of the shorter axes of the plane, logr e and (b) along another axis of the plane, 
log Mbh • 
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Fig. 2. — The curve shows the angle between (1) the best-fit mass FP, and (2) the plane 
derived using the fundamental plane in Equation @] and the relation M B h = ot log a + f3 for 
a range of values of a. The filled circles indicate the angle for actual values of a taken from 
the literature (see Tremain et al. 2002). The typical error in the measured values of a and 
the derived angle between the planes is shown at the top right in the plot. 
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log M BH (a) log M BH (L) log M BH (a) 



Fig. 3. — Comparison of black hole mass estimated using the mass fundamental plane, 
M B n(mass FP), with (a) the mass M Bii (a) estimated using Equation [1] and M B h(o") and 
(b) with the mass M Bil (L) estimated using Equation [2j In panel (c) we compare M B h(c) 
with Mbh(-^)- In each panel the dark line indicates the linear fit to the points shown, while 
the dashed line has slope unity. 



